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Introduction 



1 . 

The rising fuel costs put continuous pressure on designers to 
develop prime movers with reduced specific fuel consumption and 
smaller weight. In the case of small or medium sized gas turbines 
the demand for low weight focused increased attention on the cen- 
trifugal compressor with the requirement for improved efficiencies 
to achieve the reduction of specific fuel consumption. In the case 
of reciprocating engines the demand for reduced weight and reduced 
fuel consumption led to new concepts, like the turbo-compound and 
the adiabatic turbo-compound [1] as well as to increased interest 
in turbo-charging. These, in turn, also increased the attention 
paid to centrifugal compressors. 

Both the gas turbine and the reciprocating engine, have to be 
efficient over a wide spectrum of points of operation. If variable 
compressor geometry is to be avoided (because of increased costs, 
weight, reliability and maintenance problems) the compressor must 
have a wide range of operation in which efficiencies are good, and 
surge or choke do not occur. This demand contradicts the require- 
ment for high stage pressure ratios. 

The evaluation of the optimal compromise between maximum pres- 
sure ratio, maximum efficiency at the design point and wide range 
of operation with fixed geometry, is a difficult task. Only deep- 
er understanding resulting from improved capabilities for experi- 
mental observation and better predictive techniques for the intern- 
al flow will enable a fulfillment of this task. 
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It is the objective of this paper to review the recent develop- 
ments and the present state of the art of centrifugal compressor in- 
ternal aerodynamics. In view of the great amount of published mater- 
ial selectivity is unavoidable. Therefore the paper is focused on 
centrifugal impeller flow. Inlet and diffuser flows, as well as 
flow in axial machines, are mentioned only as much as they affect 
the centrifugal impeller flow, or are connected to its flow prob- 
1 ems . 

2 . Experimental Observations 

For a long time centrifugal impellers were designed with the 
concept of a full passage flow. The classical approach of infin- 
ite number of blades of zero thickness combined with a correcting 
slip factor is an example of this approach. However, the full pas- 
sage flow concept was already challenged as early as 1923 by Alfred 
Carrard [20] who calculated a neutral zone in the impeller passage. 

He was probably the first to propose a calculation method for this 
type of flow which about three decades later was called the jet- 
wake flow. 

The full passage concept could survive for a long period only 
because experimental insight into the impeller flow was hampered 
by the fact that it is, unfortunately, rotating. The wake in a 
rotating impeller was visualized or detected by a number of in- 
vestigators [3, 4, 5, 6] but no detailed velocity mapping was 

available until 1957, when Acosta and Bowerman [7] measured the 
flow field inside a centrifugal pump impeller with backswept 
blades (3 = 23.5°) using manometer tubes which were rotating with 
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the impeller. They did not detect a jet-wake flow, which probably 
did not occur in the type of strongly backswept impeller which they 
used. Their results, therefore, did not challenge the full passage 
flow concept, A year later, in 1958, Fujie [8] published results 
measured inside a radial exit ($2 = 90°) centrifugal compressor 
impeller using a low speed "pneumatic slip ring" , or pressure trans- 
fer cell. Fujie f s experiments were carried out at a number of tip 
speeds ranging from 37 m/sec to 105 m/sec. Although "pockets" of 
lower relative velocities were detected in the vicinity of the pas- 
sage suction side, no well developed wake region, like those found 
by later investigaters, was detected. This is surprising, but 
furthermore, Fujie’s comparison of the measurements to calculated 
inviscid irro tational and incompressible velocity profiles (Fig. 14 
of ref, 8) revealed that viscous effects now recognized and fully 
appreciated, were smaller in his results than one would have ex- 
pected. 

Ten years passed, until in 1968, Fowler [9] published results 
taken in a radial exit impeller at low speed (30 to 60 r.p.m.), 
Fowler used hot wire anemometers for his measurements. The impel- 
ler had to rotate very slowly because a technician, rotating with 
the impeller in a special cage, monitored the readings. The tip 
velocity in this experiment is estimated to have been 2.5 m/sec 
to 5 m/sec (Fowler did not publish the exact geometry and diameter 
of his impeller), but despite this low tip speed, Fowler detected 
a flow having the main features of the jet-wake structure. In fact 
he measured, in all th ro ugh- f lows from fully open flow to 50% 
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reduced flow rates, a suction side velocity that was much lower than 
pressure side velocity, - in contrast to inviscid flow predictions. 
Fowler did not detect the shear layer which separates the wake from 
the jet, and in which the velocity gradients are exceptionally 
large. However, his measurements showed that for radial exit cen- 
trifugal impellers, even at very low speeds, the inviscid flow 
models were not able to describe the real flow with sufficient ac- 
curacy . 

In 1973 J. Moore [10] published his results obtained in a radial 
exit rotating single passage. The results were reported in detail, 
including secondary flows, and were the first published data that 
describe the shear layer between a jet and a wake in a rotating 
passage. They were detailed enough to allow a quantitative analy- 
sis of a jet wake flow. The jet wake structure observed by 
J. Moore was well developed despite the small peripheral tip 
speed of about 20 m/sec and the fact that the fluid was pushed 
through the passage. Moore’s experiments were carried out at a 
number of flow rates, and it is interesting to point out that in 
contrast to that of Fowler, Moore’s jet-wake structure weakened 
considerably as the flow rate was reduced until it disappeared 
completely and wake filled up with healthy flow. Unfortunately 
this weakening of the jet-wake flow with flow rate was not dis- 
cussed in detail, as it was not consistent with Fowler’s findings 
and with the observations of later investigators. A possible 
explanation could be that Moore reduced his flow rate more than 
others did . 
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In 1975 Howard and Kittmer presented measurements carried out 
inside a backswept centrifugal impeller (6^ = 22°) with and without 
a rotating shroud [11]. Their impeller pumped water, and the veloc- 
ity field was measured with a miniature hot film probe. Like 
Acosta [9] and in contrast to Moore [10] and Fowler [9], Howard and 
Kittmer did not detect a jet-wake flow though all these experiments 
were carried out under imcompr es s ib le low speed flow conditions. 

The full passage flow in Howard’s experiments was present in both 
the closed and the open impeller versions, an indication that in 
this case the blade tip clearance flow did not affect the flow 
stability. Acosta [7] and Howard [11] measured in back swept im- 
pellers and they did not detect a jet-wake flow. Fowler [9] and 
Moore [10] measured in radial exit impellers and they found a 
jet-wake structure. Is this observation significant, taking into 
account that all these experiments were low speed and incompres- 
sible? 

Later, in 1976, Eckardt published his paper [12] on the optical 
measurement of flow in a high speed open radial exit centrifugal 
impeller, close to its design point. The tip velocity was about 
300 m/sec. Eckardt used Schodl’s L2F technique [13] and was able 
to obtain the internal impeller flow field details. His results 
showed the development of a well established jet-wake structure 
from an almost uniform inlet flow field into a highly distorted 
exit flow field. Unfortunately, the geometry of Eckardt ? s impel- 
ler was not published by him, so that his results could not be used 
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for analysis by other investigators. Luckily the missing informa- 
tion was furnished by J. Moore [14] who reconstructed Eckardt ’s 
flow passage using available data. 

Also in 1975, Mizuki et al published a very detailed report [15] 
on measurements in an open radial exit centrifugal impeller at a 
number of points of operation. The impeller rotated at a tip speed 
of 60.8 m/sec and readings were taken inside the impeller with 
4-hole yaw probes and taken out by means of mechanical seals. 
Further, static pressures were measured along the stationary shroud 
with high frequency transducers and gas velocities were measured at 
the impeller exit with a hot wire anemometer. Their results are 
interesting because like Fowler [9] and Moore [10] they extend over 
a relatively wide range of flow coefficient from y = 0.11 to 
V? = 0.33. Unfortunately, these points are not identified on the 
compressor map so that their relation to the choke and surge lines 
is not known. Like previous investigators Mizuki et al identified 
a jet-wake structure, the wake being located as usual close to the 
suction side-shroud corner, though a little displaced from the 
boundary. As in Fowler’s experiments, and even more so in Moore’s, 
the intensity of the wake was decreasing with reduced through-flow 
probably because of a general reduction of intensity of the flow 
field. This was verified by Mizuki both in his internal measure- 
ments at the impeller exit and his stationary measurements in the 
diffuser immediately after the impeller exit. An important conclu- 
sion to be drawn from Mizuki’s experiments is that the jet-wake 
flow is verified to persist over a wide range of operation. An- 
other interesting result from Mizuki’s work was that a separation 
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bubble was identified in the impeller inlet near the shroud. The 
bubble did not exist at the high flow rate (<£ = 0.33) but was 
initially detected at medium flow rates = 0.22) and became more 

intense at low flow rates (<£ = 0.11). The separation bubble was 
not detected at all by Eckardt probably because it did not exist in 
his apparatus because of the inlet geometry, or because he operated 
at a higher flow rate. 

The passages which were mapped by Fowler [9], J. Moore [10], 
Eckardt [12] and Mizuki [15] were all radial at their exits. The 
jet-wake structures detected in these passages were likely the re- 
sult of a complicated non linear combination of the effects of im- 
peller speed, exit angle of the blades, curvature of the blades, 
number of blades, existence of a shroud cover, rate of flow decel- 
eration through the impeller and flow rate. The effect of each of 
these parameters certainly can not be explored in a few experiments. 
However, a comparison of Fowler [9], Moore [10], Eckardt [12] and 
Mizuki [15] to Acosta [7] and Howard [11] suggested that the jet- 
wake flow was most dominantly affected by the sweep-back of the 
blades. This argument prompted Adler and Levy [16] to measure in- 
side a closed centrifugal impeller with swept back straight blades, 
in which a full passage flow was anticipated to occur at the design 
point. To show this was the case (and to explore the feasibility 
of optical measurements through a rotating window) they carried 
out laser doppler measurement in a closed backswept centrifugal 
impeller. They did not find a jet-wake structure. The impeller 
displayed, as anticipated, a full passage flow. This result led 
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them to hypothesize that the jet-wake flow is weakened up to total 
disappearance as both impeller speed and blade sweep-back angle are 
reduced. This view is represented graphically in Fig. (1). As 
the boundary between the full passage flow and the jet-wake flow 
is dominantly though not solely influenced by the tip speed and the 
sweepback angle a border region, rather than a border line is shown 
in Fig. 1. 

Adler and Levy's experiments were the first in which laser 
doppler measurements were carried out through a rotating window. 

They showed that a rotating window is not a severe obstacle for 
laser doppler velocimetry though special attention must be given 
to window glass quality, flatness and cleanliness. 

All the experiments described above concentrated on the mea- 
surement of the flow field inside the rotating impeller. Some 
insight into the flow inside the impeller can also be gained, 
though to a much more limited extent, from measurements at the im- 
peller exit; and in the case of open impellers, from measurements 
along the stationary shroud. Examples of this approach, which re- 
sults in a somewhat simpler experimental set-up can be found in 
the measurements of Eckardt [17] and Senoo et al [18]. Eckardt [17] 
carried out measurements at three points of operation of the stage 
investigated in [12]; one on the choke line, the second at peak ef- 
ficiency and the third about half way between peak efficiency and 
surge. He used high frequency pressure transducers on the open 
stationary shroud as well as a high frequency directionally insen- 
sitive impact tube and a hot wire directional probe after the im- 
peller exit. His measurements showed again that the jet-wake flow 
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at the impeller exit persisted over the entire operation range in- 
vestigated, though it was somewhat reduced with reducing flow rate. 
Further, a possible separation immediately behind the inlet near 
the suction side was detected. But in contrast to Mizuki [15] this 
separation bubble appeared not at low flow rates, but rather at high 
flow rates. These two results are contradictory and require addi- 
tional investigation. A third, interesting result of this work was 
the mapping of total pressure fluctuations at the impeller exit. 
Inside the wake the fluctuations reached a relative value of 14% to 
18% while in the jet they were less intense, having a value of 10% 
to 12%. This may throw some light on the distribution of losses 
and turbulent viscosity inside a rotating passage. 

A recent work was published by Senoo et al [18]. Using sta- 
tionary shroud taps to analyze the flow in a supersonic impeller, 
they were able to establish the shock wave pattern inside the im- 
peller. In addition, they discovered a new longitudinal slip band 
occurring in the subsonic flow regime. This slip band was located 
in the middle between the blades, and its thickness increased with 
flow rate reduction. It became thinner as flow rate was increased 
until it disappeared altogether. Senoo did not give an explanation 
for this flow behavior, nor did he comment on its effect on compres- 
sor performance. Evidently this newly discovered slip band should 
be further explored in more detail using optical techniques, and 
correlated with impeller characteristics. 
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The experimental observations reviewed here lead to a number 
of imported conclusions: 

1. The larger the rotor exit flow angle C^) anc * t * ie 

the rotor speed, the more likely is the development of a 
jet-wake flow towards the impeller exit. 

2. The jet-wake flow when present at the design point, exists 
also over a wide range of operation points. 

3. The jet-wake flow departs so much from inviscid flow pre- 
dictions that inviscid flow calculation are representative 
only for small 6 ^ impellers, or in case of large $ ^ 
impellers, for inlet and mid passage regions only. 

4. Small ^2 impellers generally deliver a more uniform flow 
field into the diffuser. This could explain their better 
efficiencies and wider range of operation between choke 
and surge. 
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3. 



Inviscid Flow Models 



The flow in conventional centrifugal impellers, as described in 
the previous section, is extremely complex. Clearly, the complexity 
is partly the result of not being able to design the impeller aero- 
dynamically such that the real fluid behavior is properly accounted 
for. Improved aerodynamic design procedures and better performance 
predictions, however, both require the development of techniques to 
properly analyze the flow through the impeller. Such an analysis 
must contend with the adverse pressure gradient which is imposed on 
the flow, and the presence of turbulent viscosity variations which 
are strongly influenced by mean flow conditions, curvature of the 
stream lines and the Corriolis acceleration. As will be discussed 
later, a comprehensive mathematical model describing reality is too 
difficult to solve with the current state of the art. The simplest 
approach certainly, is to calculate the flow field in the impeller 
by assuming that the fluid is inviscid. Even with this assumption 
the task is formidable, and raises the immediate query as to whether 
the effort is justified in view of the big differences obtained some- 
times between inviscid predictions and measured results [9, 10, 12, 

18]. The answer however is affirmative: inviscid solutions can be 

useful in describing backswept impeller flow where the jet-wake flow 
does not exist at the design point, in describing the upstream re- 
gions of radial exit impellers and generally to predict the inviscid 
uncorrected core of any intra-impeller flow. 

Existing methods for the prediction of inviscid centrifugal im- 
peller flows are classified into four groups I) Solutions on the 
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Hub- to-Shroud (H-S) stream surface, II) Solutions in the Blade-to- 
Blade (B-B) stream surface, III) Transonic B-B Solutions and IV) Three 
dimensional (3-D) solutions* In each of these groups models can be 
reclassified into 1) streamline curvature methods, 2) partial differ- 
ential methods . 

1) Streamline curvature methods are based on a first order 
ordinary differential equation describing the force balance in the 
direction of the quas i-or thogonals to the streamlines* The equation 
is solved along these quas i-or tho gonals marching forward along the 
passage from inlet to exit. The continuity equation is integrally 
satisfied along quas i-or thogonals . This is basically a non ellip- 
tic procedure. The elliptic character of the flow is retained some- 
how by the curve fitting of the streamlines which link upstream to 
downstream points along the passage. The curve fitting also deter- 
mines the streamline curvature. The accuracy of the results and the 
ellipticity of the solution depend strongly on which curve fitting 
method is used; this is the "non natural" or "artificial” ingredi- 
ent of the streamline curvature technique. The non-ellip ti t ity of 
the technique also lies in the indirect manner in which downstream 
boundary conditions and backpressure are satisfied, and in the way 
their upstream effects on the entire flow field are predicted. 

An early and extremely simple example of this approach is the 
method of Adler and Ilberg [19] where the assumption of linear 
variation of streamline curvature from hub to shroud along the 
quasi-or thogonals allowed integration of the equilibrium equation 
in the direction of the quas i-or thogonals . The method enabled the 
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flow field to be obtained without the use of a computer. 

2) Partial differential methods are based on the solu- 
tion of the partial differential equation of the inviscid steady 
flow. Momentum equilibrium and continuity are both satisfied at a 
point, and the elliptic nature of the subsonic flow is well pre- 
served both in the form of the partial differential equation and, 
consequently in the way the boundary conditions can be specified. 

The formulation is based in most cases on the classical work of 
Wu [20] or on equations similar to those given by Wu. Wu defined 
"special derivatives" valid only on H-S or B-B stream surfaces. 

These special derivatives allow an enormous simplification of the 
problem by splitting the 3-D formulation into two mathematically 
"two dimensional" formulations. These are, however, interconnected 
by the "stream sheet thickness" thus retaining the three dimension- 
ality of the physical situation. The partial differential equations 
can be solved either by finite differences or by finite element 
methods. These will be described later. 

After having briefly explained the streamline curvature and the 
partial differential approaches let us continue the review along 
the classification into H-S, B-B and 3-D groups. Due to the large 
number of methods known, only a brief description can be given here. 
After this the various methods will be compared to each other and 
to experiments . 

3.1 H-S solutions : Historically, streamline curvature methods 

were the first to be developed. They required more modest computer 
storage and simpler computational techniques were needed because 
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they involved only first order, ordinary differential equations. 

A frequently used streamline curvature program is Katsanis ’ early 
development [21]. Katsanis expressed the equation of motion of 
the H-S flow in terms of curved coordinates based on the stream- 
lines and quasi-or thogonals . The resulting equation was solved 
simultaneously with the continuity equation in its integral form 
along the q ua s i-o r t ho go na Is . The numerical solution iterated on 
the constant mass flux between neighboring streamlines. The 
curvature of the streamlines was modified from iteration to itera- 
tion until convergence was achieved. 

An additional method based on the streamline curvature tech- 
nique was developed by Novak [22, 22]. As in the previous method 

the equation of motion was expressed in terms of the streamline 
curvature and the solution iterated on the constant mass flux be- 
tween two successive streamlines. Unlike Katsanis, Novak modified 
his initially assumed streamlines along radii, rather than along 
q uas i-o r tho gonals . This simplified the method but restricts its 
application to axial flow machines only. Davis [24] compared both 
Katsanis 1 and Novak’s methods to a method developed by himself. 

He concluded that all three versions of the streamline curvature 
method were essentially identical. A general analysis of the stream- 
line curvature approach can be found in a paper by Smith [25]. 

As mentioned above the partial differential methods are based 
in most cases on Wu’s formulation [20] or its derivatives. An 
early finite difference solution was given by Wu [26], and later 
Marsh [27] used a modified version of Wu’s model to obtain a 
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solution of the flow field in the H-S surface. Later, in 19 73 , 
Davis published his technique [28] which was essentially similar 
to Marsh’s. It was developed for axial machines but could be ap- 
plied also to centrifugal compressors in an analogous way to 
Katsanis ' and McNally’s evolutionary work. Katsanis and McNally 
developed a finite difference method [29], which in its early 
version was limited to flows up to 45° from axis. Their later 
version [30] extended the application to radial flows as well. 

Both methods, those of Davis and Katsanis and McNally, were based 
on the use of the stream function equations. The equations were 
solved with a finite difference technique, and the basic method 
was limited to subsonic flows. Locally supersonic flows could 
only be handled by a combination of the stream function solution 
and a velocity gradient (streamline curvature) technique [42]. 

A second way to solve the partial differential equations of 
the H-S flow is the finite element technique. The basic equations 
solved are essentially identical to those used in the finite dif- 
ference solutions based on Wu’s model, except for modifications 
introduced for the convenience of the technique. The first finite 
element solution was completed in 1974 by Adler and Krimerman 
[31]. This was followed in 1975 by Oates et al [32] and in 1976 
by Hirsch and Warzee [33]. A basic difference between the finite 
element and the finite difference techniques is that the discre- 
tization of the continuous flow field in finite difference form is 
geometrically restricted by the coordinate system chosen for the 
problem (unless complicated and time consuming transformations are 
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applied) while in finite element form no restriction is imposed 
on the location of the nodal points in the grid. The freedom in 
location of grid points enabled Adler and Krimerman [31] to locate 
their grid points on streamlines. This allowed a simpler satisfactioi 
of the energy equation which, along streamlines, could be used in 
its integrated form. Further, the fact that grid points were 
located on streamlines made the incorporation of the solution into 
a fully 3-D solution easier. It was with this recognition that Adler 
and Krimerman choose the finite element technique for their solu- 
tion. 

Oates et al [32] applied a variational principle, producing, 
in effect, the integral of the H-S momentum, to create a system 
of non-linear algebraic equations. The system of equations was 
solved using the Newton-Raphson technique. To formulate their 
method Oates et al used the well-known actuator disc method. The 
technique was fundamentally different from that of Adler and 
Krimerman who did not use a variational principle specially deriv- 
ed for their principal equation. Rather, they used the variational 
principle of Poisson’s equation, giving their governing equation an 
artificial shape of a Poisson equation. This in turn required an 
iterative procedure. Further Adler and Krimerman did not use the 
actuator disc approach at all. 

Hirsch and Warzee [33] used a weighted residual Galerkin pro- 
cedure. Like Adler and Krimerman they gave their governing equa- 
tion the shape of a Poisson equation, avoiding in this way the 
problem of deriving a variational principle for their governing 
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equation. They also introduced NASA correlations into their tech- 
nique to specify losses and deviation angles. It should be pointed 
out here that all these techniques were initially applied to axial 
flow compressors. There is, however, no special reason to believe 
that these techniques will fail when used to predict the flow in 
centrifugal impellers. In fact, Adler and Krimerman used their 
technique [31] successfully to analyze Eckardt f s impeller [12] with 
out any difficulties. This will be discussed later. 

No systematic comparison of the simultaneous application of all 
these techniques has been reported. Further, no systematic compari- 
son of the predictions of the various methods with the results of 
experiments has appeared. There are, however, isolated comparisons 
between two or three techniques or between a numerical solution and 
an experiment to be found in literature. Davis compared his stream 
line curvature method to those of Katsanis and Novak and found them 
to be essentially identical [28]. In 1972 Frost [34] presented a 
streamline curvature technique and compared its results to Marsh’s 
method [29]. The agreement between the two techniques was very 
good but both do not agree too well with experimental results 
(see Fig. 2). Differences were observed not only in velocity mag- 
nitude and direction but also in velocity profile shapes, indicat- 
ing a fundamental problem; perhaps the effect of viscosity? 

Adler and Krimerman [31] compared their results calculated 
usinga finite element method to Wu's results [26] evaluated with 
a finite difference method. The two results agreed fairly well 
and showed the same trends in the distributions of velocity and 
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static pressure along the casing and the hub. When compared to 
experimental data both analytical techniques gave static pressure 
distributions showing trends similar to the measured profiles 
(see Fig. 3). Adler and Krimerman's results were closer to the 
measured profile when blade thickness effects and a total loss 
coefficient were included in the calculations. 

In 1974 Davis and Millar [35] reported a comparison between 
finite difference and streamline curvature methods. Their com- 
parison is illustrated in Fig. 4 where test results of a NASA 
single stage transonic compressor rotor are also shown. The two 
calculations predict similar velocity, pressure and total tempera- 
ture distributions and agreement with experiments is reasonable. 
Finally in 1976 Hirsch and Warzee’s finite element method was 
compared to experimental data obtained in a single stage NASA 
compressor (see Fig. 5). Agreement seems to be exceptionally 
good for an inviscid theory. 

A final word of explanation and caution should be added here. 
Most of the H-S flow prediction methods mentioned so far and the 
comparisons discussed were carried out for axial flow machines 
(Fig. 2-5). The methods thems elves, however , are of interest to 
centrifugal impeller flow studies (see footnote*). But it should 
always be kept in mind, that three dimensional and viscous effects 
in centrifugal impellers are considerably more significant than in 
axial rotors, especially if a jet-wake flow is present. Therefore, 



*For example the method of Ref. [31] was used unaltered in a 
centrifugal radial exit impeller. 
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in the case of conventional centrifugal impellers inviscid H-S 
techniques will give results which are further from reality than 
the comparisons in Figures 2 to 5 would suggest. Further, all 
the H-S methods reviewed so far assume axisymmetric flow. In the 
higher loaded centrifugal impellers this assumption could lead to 
big discrepancies as will be shown later in this paper. 

3.2 B-B Solutions : As in the case of the H-S surface solu- 
tion the formulation of the B-B stream surface flow solution is 
somewhat artificial in that the assumed surface does not exist 
physically in the centrifugal impeller. The B-B solutions, however, 
serve to analyze trends required to be known for passage optimiza- 
tion rather than predict exact flow fields. Further, they can serve 
as essential "building blocks” in the development of more complete 
3-D and viscous methods. Ba s i ca 11 y , t h e two groups of solutions 
on the H-S and the B-B stream surfaces have much in common. However, 
the most pronounced difference between the two is the lack of well 
defined boundaries at the inlet to, and the exit from the blade 
passages on the B-B surface. These boundaries are given on the 
H-S surface by the hub and shroud extensions in the upstream and 
downstream directions, and are therefore considered to be input 
data. In constrast, the boundaries, in form of the leading - 
and trailing edge stagnation streamlines, are determined on the 
B-B surface in the course of the solution. 

An approach to the B-B problem not used on the H-S surface 
is the method of singularities applied in all cases on a surface 
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of revolution [36, 37, 38, 39, 40, 41]. The surface is the center 

of a stream sheet which, in reality, is of varying thickness be- 
cause of continuity with compressibility. The method has problems 
at the leading edge and the trailing edge, and the ways to overcome 
them are not always general enough. Martensen’s [36] and Ogawa's 
[40, 41] methods, for example, develop difficulties when the blades 

are very thin (a common case in centrifugal impellers) or when the 
trailing edge is thick and rounded (again common in centrifugal 
impellers). It is reported that Wilkinson [37] was able to solve 
this problem; however, for an axial machine. In a centrifugal 
compressor, where viscous effects on deviation angle (and slip) 
are more significant, this problem is not yet solved. A severe 
problem in the application of the singularities technique in centri- 
fugal impellers is the strongly varying streamsheet width. Ogawa 
and Murata [41] solved this problem by approximating the width 
variation by an analytical expression and, in a more exact manner 
by a numerical solution of their equations. From their work, it 
appeared that in spite of a large stream sheet variation, the solu- 
tion yielded acceptable results. 

As on the H-S surface, the streamline curvature method can 
be applied to the B-B surface flow prediction. It has been used 
less often however, probably because of the stagnation streamlines 
problem. One of the first methods was developed by Katsanis [42]. 
The method was based on an equation for the velocity gradient 
along quasi-or thogonals between the blades and was similar to the 
method of quas i-or thogonals on the H-S surface. The solution was 
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carried out on a surface of revolution generated by rotating a 
previously calculated H-S streamline around the axis* Later, 
Katsanis used this streamline curvature technique in the supersonic 
regions of a transonic finite difference method which he developed. 
Wilkinson also developed a streamline curvature method [43]. He 
gave much attention to the location and shape of leading- and 
trailing edge stagnation streamlines, which is essential in this 
type of calcualtion. Wilkinson, who carried out basic work on 
streamline curvature techniques, their accuracy and convergence 
characteristics as affected by the curve fit procedures [44], was 
able to optimize the techniques as to the number of iterations 
required until convergence by deriving an optimum damping factor. 
Comparison of his results with experiments is given in Fig. 6. 

Most frequently used in the B-B flow problem solution are 
finite difference techniques. Stanitz, a pioneer in this field, 
published a finite difference method as early as 1948 [45, 46, 

47]. He determined the flow field by solving stream function 
equations. Stanitz also calculated off-design B-B flow fields 
and predicted pressure side separations (see Fig 7) which, as is 
now well known, do not exist because of the important effect of 
viscosity on flow field. Stanitz’s results illustrate clearly 
that inviscid methods are not able to predict, even approximately, 
off-design flow fields on the B-B stream surface. They should 
be used only at design point, and even here they will not be able 
to predict the jet-wake flow which has been measured by many 
inves t i ga to rs . 
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In 1968 Katsanis published a finite differences method to 
predict the B-B flow field [48]. Like most B-B flow calculation 
techniques, Katsanis 1 also suffered from the inaccuracy introduced 
by the sensivity of the results to the location and shape of the 
leading edge stagnation streamlines. To improve the accuracy of 
the prediction and to enable more detailed information in this 
critical part of the flow field, Katsanis and McNally [49] developed 
a method to predict the details of the flow field near the leading 
edge of the blades. They further developed their technique to 
handle splitter blades in centrifugal compressors [50]. The pro- 
grams of Katsanis and McNally are popular and widely used but 
they also have the deficiency from which all inviscid B-B programs 
suffer; they can not predict the jet-wake flow, which is essential- 
ly a viscous phenomenon. This is a severe drawback of all inviscid 
B-B me tho ds . 

To complete the list of finite difference methods, the work 
of Smith and Frost should be noted [51]. Smith and Frost developed 
a finite difference technique based on Wu’s formulation [20]. 

They programmed a solution using Marsh's scheme [27], and also 
programmed a streamlind curvature technique similar in its general 
features to those already mentioned. Using the two computer pro- 
grams, Smith and Frost compared the two techniques by carrying 
out calculations for identical problems. In the case of an axial 
compressor a fair degree of agreement was achieved, except on the 
suction side leading edge where the velocity distribution computed 
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by streamline curvature was oscilating. Smith and Frost believed 
that a finer mesh would have eliminated this problem. In the case 
of an axial turbine the streamline curvature technique behaved badly, 
and velocity oscillations were encountered both at the leading and 
the traling edges, caused probably by difficulties in evaluation of 
streamline curvatures. This problem did not exist in the finite 
difference solution. Generally it can be said that the streamline 
curvature technique required less computer storage than the finite 
difference technique, but the computing time was longer. In contrast 
to the finite difference method, the streamline curvature technique 
was not restricted by the local Mach number and could therefore be 
used in cases with local supersonic regions. 

A general problem in numerical techniques is how to optimize 
convergence, accuracy and computation time. To demonstrate that the 
convergence and execution time of early programs can probably be 
improved in many cases, one can cite Deshpande f s work [52]. 

Deshpande developed a new successive overrelaxation algorithm to 
improve Katsanis 1 original algorithm [48]. The new algorithm re- 
duced execution time by a factor of 2.5 for a typical blade config- 
uration. Further, Deshpande’s algorithm enabled convergence for 
cases in which Katsanis 1 algorithm did not converge at all. Ap- 
plication of Deshpande's new algorithm in Katsanis’ programs [48, 

49, 50] required only minor modifications and is reported by Deshpande 
to have reduced computing costs considerably. 

Davis and Millar [51] showed very clearly how complex the finite 
difference stencils can become in cases of unusual boundary geometry 
or with irregular grids imposed by the coordinate system. This 
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problem can be removed by using the finite elements method. 

B lad e- t o-b lad e finite element methods were published by Adler 
and Krimerman in 1977 [54], followed by Prince in 1978 [55]. Adler 

and Kriraerman’s approach was very similar to that taken in their 
H-S solution [31]. The main difference was that special care was 
taken to predict correctly the leading-and trailing edge stagnation 
streamlines. This was done by imposing the conditions of period- 
icity and constant angular momentum on the approaching flow field. 

The downstream stagnation streamline was evaluated using the Kutta- 
Jukowski condition and flow field periodicity. As in the H-S 
finite element solution, the B-B solution took advantage of the 
freedom to choose the location of grid points. Consequently, grid 
points were always located on streamlines, thus allowing again the 
use of the integrated energy equation. The model solved was Wu's 
B-B formulation [20], which was again given the artificial form of 
a Poisson equation to overcome the lack of a variational principle 
to Wu's principal equation. A comparison of Adler and Krimerman’ s 
results with experiments and with calculations by Stanitz, is shown 
in Fig . 8 . 

Prince [55] followed a somewhat different approach. He used a 
Galerkin variational analysis to satisfy the continuity equation for 
steady potential flow. He then used a Newton-Raphson method to 
solve the resulting system of nonlinear algebraic equations. Prince 
used continuity conditions to evaluate the upstream and downstream 
stagnation streamlines. To do this he used a distribution of stream 
sheet thickness which was required to be known a priori. Prince’s 
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calculated results compared well with experiments as is shown in 
Fig. 9. 

3.3 Transonic B-B Flows : The case of transonic B-B flow is 

treated separately because the nature of this type of flow changes 
from elliptic when the Mach number is less than unity to hyperbolic 
when the Mach number is greater than unity. The difficulty in 
transonic solutions is that elliptic flow equations are solved with 
numerical techniques which differ considerably from the numerical 
techniques used in the solution of hyperbolic partial differential 
equations. The problem becomes even more complicated because the 
boundaries between the elliptic - and the hyperbolic regions of 
the flow (on which the equations are parabolic) is not known 
a priori but must be determined during the computation. 

An early transonic solution was again derived by Katsanis [56J. 
Katsanis ’ method can be .applied to any impeller geometry (axial or 
centrifugal). The solution was obtained as a combination of a 
finite difference technique, a stream function solution and a veloc- 
ity gradient method. The finite difference solution, carried out 
at a reduced flow rate, provided information required for the veloc- 
ity gradient method. For flow which was totally subsonic the finite 
difference method was used. In supersonic parts of the flow field, 
the streamline curvature technique was employed. Katsanis used 
different numerical techniques to obtain the solution in the dif- 
ferent regions which required automatic detection of the region type 
and, accordingly, a switch over from technique to technique. These 
complications could be overcome if a single solution could handle 
both sub- and supersonic regions of the transonic flow field. 
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Although the equations for steady flow are hyperbolic for super- 
sonic flows and elliptic for subsonic flows, for unsteady flow the 
equations are hyperbolic both for super- and for subsonic flows. 
Therefore, if the time dependent equations of continuity, momentum 
and energy are used, a single solution technique is applicable for 
both types of flow. The steady flow, including both super- and sub- 
sonic regions, can then be regarded as the final steady state of 
the time dependent flow. This concept leads to a unified method used 
identically in the two flow regimes. Generally an approximate flow 
field is used as an initial condition and then the equations are 
integrated forward in time until the final steady solution is reached 
The technique is known as the "time marching” method. In some cases, 
severe stability problems are encountered. They can often be solved, 
however, by reducing the time steps or by using an artificial vis- 
cosity which has a damping, effect. 

In 1971, Marsh and Merryweather [57] described a stable time 
marching technique which was based on finite differences and did not 
rely on the use of artificial viscosity to achieve stability. It 
was found that several stable procedures could be developed for 
flow in convergent-divergent nozzles. The characteristic feature 
of the stable schemes was that the derivatives of all quantities 
other than pressure were approximated by backward differences, 
while the derivative of pressure contained a forward element. The 
computer program developed by Marsh and Merryweather [57] was rela- 
tively slow, some 2700 iterations being required to obtain a solu- 
tion to an accuracy of 0.01 per cent. Further work by Daneshyar 
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and Glynn [58] was based on the method of characteristics which led 
to a much faster method of calculation. The method has been ex- 
tended by Glynn to deal with cascade flows. 

In 1971, McDonald [59] used a time marching method to calculate 
the pressure distribution around aerofoils in cascade. The problem 
was formulated in terms of a finite area approach which led to the 
conservation equations in an integral form. The flow was assumed 
to be isentropic on the grounds that only weak shocks were normal- 
ly encountered in cascades. McDonald obtained very good agreement 
between his calculated pressure distribution and that measured ex- 
perimentally in a cascade. The use of the isentropic flow assump- 
tion is interesting in that Marsh and Merryweather had tried the 
same assumption for purely subsonic flows and had experienced a 
severe numerical instability. In impellers, where the flow is cer- 
tainly not isentropic, these methods [57, 58, 59] have to be modi- 

fied. 

Also in 1971, Gopalakr ishnan and Bozzola [60] published a time 
marching transonic technique. They used a predictor-corrector, 
two step, time splitting method to solve the partial differential 
formulation of the transonic flow problems. The method had sta- 
bility problems unless the time step was smaller than a value pre- 
dicted through a linear analysis, but even this was not always 
sufficient to give stable computation. To overcome this difficulty, 
a spatial smoothing procedure could be used. Excessive smoothing, 
however, caused inaccuracies in the results, as would be expected, 
and the problem then is the prediction of the negative effect of 
the smoothing. 
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In 1974, Denton [61] proposed a time marching scheme for cas- 
cade flows using a simpler grid than that of McDonald. Denton’s 
grid consisted of qua s i- s t r e ami ine s and straight lines across the 
blade passage. The conservation equations for mass, momentum and 
energy were derived for a control volume. Instead of assuming 
isentropic flow, Denton assumed constant stagnation enthalpy, 
an assumption which became exact when the solution converged to the 
steady state flow. In Denton’s scheme, the pressure at the central 
point of an element was assumed to act on the upstream face of the 
element, whereas the velocity at the center controlled the flow 
through the downstream face. The maximum time step for this scheme 
was found to be larger than for the method of Marsh and Merryweather 
[57] . 

Denton has applied his time marching method to calculating the 
blade- to-blade flow in several cascades and has obtained encouraging 
results. He has also extended the method to three-dimensional flows, 
although this did require a large amount of storage in the computer. 
The predictions obtained with this program have been compared with 
experiments performed with a rectangular duct having 60° of turning. 
Good agreement was obtained between the calculated and experimental 
pressure variations for the four corners of the duct. It should be 
possible to extend the time marching scheme to deal with three- 
dimensional flow in cascades. 

It is noted again that Denton’s technique and some of the other 
methods mentioned were developed for axial flow stator cascades. 

The methods can, in principle, be modified to handle the flow in 
centrifugal impellers. It simply has not yet been done. A 
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computational technique to solve the flow in axial rotors was devel- 
oped in 1975 by Kurzrock and Novick [62] in an approach similar to 
that of Gopalakr ishnan and Bozzola [60]. They solved the Navier- 
Stokes equations using the time marching technique. The method was 
based on a time splitting finite difference technique, using a 
two step numerical evaluation of the time derivatives of the partial 
differential equations, together with a predictor-corrector tech- 
nique. The spatial derivatives were numerically evaluated by back- 
ward differences for the predictor step and by forward differences 
for the corrector step. Figure 10 illustrates the good agreement 
found between Kurzrock and Novick' s results, and experimental 
measurements . 

A completely different approach to the problem was introduced 
by Sobieczky, Fung and Seebass L 6 3 ] . Their method was developed 
for single airfoils but can be extended to stationary or rotating 
cascades [64]. Sobieczky et al developed a technique which enabled 
the design of shock free transonic aerofoils, in principle elimin- 
ating the losses associated with shocks and shock-boundary layer 
interactions. Their procedure was based on the use of any numerical 
method for subsonic, compressible isentropic flow. The algorithm 
of the selected method was modified so that if the flow became hyper- 
bolic, the basic equations in that region were altered, using a 
"fictitious gas" concept, such that the system of equations reverted 
to elliptic behaviour. In this way the complete flow field over a 
given configuration was first calculated using an elliptic procedure. 
The calculation served to define sonic surfaces. 
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Outside the sonic surface the solution obtained satisfied the correct 
unmodified equations, and the potential at infinity had the correct 
value for the circulation. The second step in the procedure was 
to modify the surface geometry inside the sonic surface, using a 
hyperbolic procedure for the real gas, to generate a shock-free 
supersonic flow with the same conditions along the sonic line. 

Here a problem arises as it is well known that shock free, two or 
three dimensional irrotational near-sonic flows are physically and 
mathematically isolated. Any small change in flow - or boundary 
conditions might in practice lead to the formation of shock waves. 
However, if the shock waves are very weak the method will have prac- 
tical significance. While Sobieczky ! s method was developed for 
inviscid flow, the same procedure can be employed iteratively with 
a boundary layer calculation to achieve shock free viscous designs. 

3.4 3-D Inviscid Solutions : The prediction techniques pre- 

viously described were all two dimensional. Their relevance to 
describing real impeller flows must be questioned, especially in 
the case of highly loaded centrifugal impellers where three dimen- 
sional effects are mush more significant than in axial rotors or 
cascades. The answer to this question can be obtained in two ways: 

1) comparison to experiments, 2) comparison between 2-D and 3-D 
solutions both carried out for a given impeller at the same oper- 
ating conditions. First, currently available methods which can 
predict 3-D centrifugal impeller flow (or related axial rotor flow) 
will be reviewed. 
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An early attempt to predict the 3-D inviscid impeller flow was 
made by Schilhansl [65]. The work was restricted to incompressible 
flow and was developed with the severely restrictive assumption that 
all stream surfaces were surfaces of revolution. The method was 
based on cascade theories then available in 1965. In fact the method 
of singularities used by Schilhansl restricted the application even 
further to irrotational flow, a severe limitation in turbomachinery 
flows. No results were reported in Schilhansl's paper [65], but 
the work illustrated very clearly what the difficulties were in 
the solution of 3-D flows. Six years later, in 1971, Katsanis pub- 
lished a 3-D compressible flow prediction method [66J. Katsanis 
used the ordinary differential equations of the streamline curvature 
technique, which define the values of velocity gradients along quasi- 
orthogonals. The continuity equation was integrally satisfied as 
is normally done in the streamline curvature methods. Katsanis 
solved the H-S and the B-B velocity gradient equations simultane- 
ously, with the condition that either the weight flow was specified 
or that the flow was choked, to determine the velocity distribution 
over the blade surfaces. The method incorporated a number of assump- 
tions, some of which could affect considerably the accuracy of the 
results in extreme geometries or flow conditions. The most re- 
strictive assumptions were, 1) that there was either a linear vari- 
ation of curvature or of radius of curvature along an orthogonal, 

2) that there was no change in radius along a B-B orthogonal, and 

3) that the flow angle along a B-B orthogonal was constant. 

In 1972 Senno and Nakase published a quasi 3-D method [67]. 
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The method included again the limitation that the b 1 a d e - t o -b 1 a d e 
stream surfaces were a xi s ymrae tr i c . (In reality these stream sur- 
faces are highly twisted and distorted.) Senoo and Nakase’s method 
was based on two interacting procedures; one, a streamline curva- 
ture H-S procedure and the other a partial differential blade-to- 
blade technique developed by Senoo and Nakase a year earlier [68]. 

To achieve a final solution the method iterated from H-S to B-B 
and back, but the B-B stream surfaces were unaltered during the 
process and remained surfaces of revolution. The H-S method was 
rather conventional while the B-B solution was based on Prasil 1 s 
independent variable transformation. The transformation simplified 
the method, but required the B-B stream surfaces to be axisymme t r ic . 
A comparison between Senoo and Nakase’s calculated results and 
experiments is given in Fig. 11. 

Six years later, two further works were published; Bosnian and 
El-Shaarawi [69], and Novak and Hearsey [70]. Bosman and El- 
Shaarawi used Wu T s approach to provide the mathematical model. 

They worked iteratively with the solutions obtained alternately 
on the H-S and B-B stream surfaces. However, the B-B surface was 
again always a surface of revolution and the H-S surface was cal- 
culated for an averaged mass flux. In reality the H-S surfaces 
and the mass flux through them change considerably across the 
passages. Thus their method was not truly three-dimensional. The 
solution of the principal equations was obtained using Marsh's 
matrix inversion method rather than by a relaxation method, because 
comparative studies carried out by Hill [71] and Bosman [72] showed 
that the matrix inversion method was considerably faster than 
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the relaxation method. No comparison with experimental results 
was given in Bosman and E 1-Shaar awi * s paper. But the differences 
between a 2-D computed flow field and their quasi 3-D results 
are demonstrated to be significant enough in the case of a typical 
centrifugal compressor impeller to justify the more elaborate quasi 
3-D computation. 

Novak and Hearsey [70] also used two 2-D methods, one on the 
H-S surface and the other on the B-B surface to generate quasi 
3-D results. Both 2-D component methods used the streamline curva- 
ture technique. The two programs were coordinated to yield a quasi 
3-D solution, but as in the previous example, the B-B surfaces were 
surfaces of revolution, and the H-S surfaces were stream sheets 
roughly parallel to the blades. In other words, the H-S surfaces, 
which in reality vary from pressure side to suction side of the 
rotating passage, are represented in Novak’s and Hearsey f s work by 
a single surface. Comparisons with experimental measurements were 
given for an axial nozzle passage and an axial compressor cascade, 
and agreement was good. Centrifugal compressor flows were not an- 
alyzed, and therefore the accuracy of the technique in this case of 
more distorted flow is not known. 

In 1978 two new 3-D inviscid methods were published, Hirsch and 
Warzee [73] and Krimerman and Adler [74] developed quasi 3-D methods 
based on Wu’s approach. The method of Hirsch and Warzee [73], again 
used H-S and B-B formulations which were combined iteratively to 
yield a quasi 3-D flow. In Hirsch and Warzee’s method the interac- 
tion between the H-S and the B-B solutions was stronger than in the 
previous methods because the stream sheet thickness and the streamline 
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angles were transferee! from the H-S solution to the B-B solution. 
Further, the fluctuation terms (as defined by Hirsch and Warzee [33]) 
were determined by a number of successive B-B solutions and then 
transfered to the H-S calculation. This was a step toward a more 
three-dimensional calculation; but nevertheless the B-B stream 
surfaces were always axisymmetric and the H-S solution, although 
corrected by computed fluctuation terms resulting from the B-B 
results, was still carried out on a single surface. As in their 
first method [33] Hirsch and Warzee used the finite elements tech- 
nique for their solution. 

Krimerman and Adler f s method [74] was based on H-S and B-B 
solutions which they developed earlier [31, 54], and used the finite 

elements method of solution. However, in contrast to the previous 
methods in Krimerman and Adler’s approach the B-B surface was not 
required to be axisymmetric and the real stream sheet thickness was 
taken into account as a spatially three-dimensional function. 
Further, the H-S flow was not represented any more by a single 
stream surface but was calculated on an arbitrary number of H-S 
surfaces which varied in shape from pressure side to suction side 
across the passage. Krimerman and Adler solved Wu 1 s model, which 
they modified slightly. They began by calculating a single mean 
H-S solution. Then they calculated a number of axisymmetric B-B 
solutions between the hub and the shroud. As a third step, a number 
of H-S solutions were carried out, giving the circumferential vari- 
ation of the flow field from pressure-to suction side. The fourth 
step was to calculate a number of B-B solutions in which the real 
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variation of the H-S flows from pressure side to suction side of 
the passage was accounted for. Steps three and four were itera- 
tively repeated until convergence was achieved. In this method 
[74] all the quantities and stream surface geometries were truly 
three-dimensional. No restricting assumptions about stream surface 
geometries are made and no pitch averages or fluctuation terms 
were required. The results were restricted only by the inviscid 
nature of the fluid and by the fact that corner streamlines were 
fixed by the boundary conditions to be identical with the geo- 
metrical passage corners. In reality streamlines can wrap around 
the geometrical corners of the passage. Figure 12 shows a sample 
of the results obtained by Krimerman and Adler for the flow in a 
centrifugal impeller. In Fig. 12, the intersection lines of quasi- 
orthogonal planes and the computed H-S and B-B surfaces are shown 
for four selected planes. The departure of the B-B surfaces, 
computed by the 3-D method from the axisymmetric B-B surfaces 
used previously is seen to be considerable. Figure 12 also shows 
the variation of the H-S stream surfaces from pressure side to 
the suction side of the passage. In previous methods this varia- 
tion was either not taken into account or it was represented by 
pitch averages and fluctuation terms. Figure 13 shows how the 
method converged from the initial steps to the final 3-D solution. 

The comparisons between the predictions of various inviscid tech- 
niques and the results of Figs. 3, 4, 5, 6, 8, 9, and 10 are encour- 
aging. Only in Fig. 2 and Fig. 11 are discrepancies clearly evi- 
dent. The comparison shown in Figure 14 fits well into this en- 
couraging scene except for the data near the impeller exit, where 
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the measured average pressure is far from the prediction. Is this 
a failure of the analysis or a measurement error? This question 
must be examined with greater care. 

Adler and Kriraerman [75] also used the technique described in 
ref. [74] to calculate the flow in Eckardt ’ s impeller [12], using 
the geometry reconstructed by Moore [14]. The results for some 
selected blade- to-b lade sections in the center of the pasage are 
given in Fig. 15. The predictions appear to be good except perhaps 
near to the exit. However, Fig. 16 gives the comparison between cal 
culations and measurement in the same impeller but close to the 
shroud where Eckardt detected the wake. Here the prediction is 
seen to be poor. The reason for the poor agreement between calcu- 
lation and experiment seen in Fig. 16 and at the exit regions in 
Fig. 14 and Fig. 15, is almost certainly the effect of viscosity 
which was eliminated in the analytical models. 

The failure of the inviscid models to predict the flow in 
centrifugal impellers is manifested when a jet-wake flow is present. 
If the classification of Fig. 1 is supported by additional experi- 
ments it could be used also to specify classes of flows for which 
inviscid models are acceptable and those for which inviscid models 
fail. The situation here is quite different from that found in 
axial flow machines and is reminiscent of the inadequacy which is 
found in applying inviscid models in cases of separated flows. The 
wake is probably similar to a separated flow, although different in 
that it must be influenced by the strong centrifugal and the Corioli 
acceleration fields which are present -in rotating passages. These 
features will be discussed in the next paragraph. 
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Another major limitation of the inviscid models described in the 
previous paragraphs is their inability to predict correctly second- 
ary flows, which are to a large extent dominated and initiated by 
viscous effects. The secondary flows, so well described by Eckardt 
[12], can have a significant effect on the impeller flow field and 
consequently on the impeller performance. Secondary flow effects 
must be included in a really representative prediction method. 
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4 . The Structure of Turbulence and the Turbulent Viscosity in 

Rotating Passages 

A conclusion in the previous paragraph was that viscosity can 
be a dominant influence in centrifugal impellers and that in many 
cases inviscid methods must fail to predict reality. The role of 
viscosity is influential, not only because of the adverse pressure 
gradients present but also because viscosity variations caused by 
curvature and Coriolis accelerations can be large. Because the 
flow in most practical impellers is turbulent, work on laminar 
viscous flows in centrifugal impellers are not in general reviewed 
here. One result however is noted. Grundmann’s study [76] showed 
that in laminar flow, centrifugal and Coriolis acceleration terms 
and the additional curvature terms had a strong influence on the 
location predicted for the separation point using boundary layer 
theory. If the centrifugal acceleration exceeded a certain limit, 
the solution scheme failed. This makes the use of boundary layer 
theory doubtful for this class of flows. 

Most practical impeller flows are turbulent, and the pattern of 
these flows is, among other factors, the result of significant tur- 
bulent viscosity gradients. The aim of this paragraph is to review 
the main reasons for the viscosity gradients, the nature of vis- 
cosity in rotating passages and the effect of viscosity gradients 
on the flow field. The literature on this subject has increased 
considerably in recent years. On the basis of reported studies 
it appears that the viscosity variation in the flow field inside 
a centrifugal impeller is the result of the interaction between 
a body force field and the velocity field. The viscosity variation 
depends on the relation between these two fields and is greatest 
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at locations where velocity gradients are large. For the purpose of 
this review the literature will be classified into three groups: 

1) effects of curvature, 2) effects of rotation, and 3) combined 
effects and application. 

4.1 Curvature Effect : In 1969 Bradshaw described an analogy 

between streamline curva t ur e-and buoyancy effects on turbulent 
viscosity [77]. He showed that there was a similarity between 
thermally stratified turbulent viscous flows, (which were by then 
better explored, especially in meteorological problems) and turbu- 
lent viscous flow along curved streamlines. He also discussed, 
briefly, effects of rotation. His major conclusion was that the 
effects of curvature, either concave or convex, on the viscosity 
were large if the thickness of the shear layer exceeded l/300th 
of the radius of curvature of the streamline. These effects in- 
creased with increasing Mach number. Later, in 1973 and 1975, 

So and Mellor studied experimentally the effects of curvature on 
turbulent viscosity [78, 79]. They found that, along a convex 

wall, the Reynolds stress was decreased so strongly that it van- 
ished in the middle of the boundary layer. In other words, they 
found that the boundary layer was ’’laminarized 11 and that turbu- 
lence was suppressed. In contrast, they found that along a concave 
wall the turbulent intensities inside the boundary layer were 
substantially increased. These phenomena have a striking effect on 
centrifugal impeller flows. In 1975, So [80], derived a formula 
for the variation of the turbulent velocity scale with the Richardson 
number for curved shear flows. He used the Reynolds stress equations 
and assumed that production of turbulent energy balanced the viscous 
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dissipation. This provided a first quantitative account of the 
phenomenon. In the same year Irwin and Smith [81] also published 
a quantitative prediction of curvature effects. They verified 
that even small amounts of streamline curvature have a surprisingly 
strong effect on the eddy viscosity. Using data from curved 
boundary layers, they were able to predict the curvature effect 
properly if the curvature terms were included in the model. It 
is important to note the considerable effect of even small curva- 
tures. It can explain some observed impeller flow field character- 
istics, such as suction side separation and pressure side attach- 
ment, which are contrary to inviscid theory predictions. 

Recently Shivaprasad and Ramaprian published two experimental 
studies [82, 83] which again verified that turbulent viscosity was 

considerably affected by an even small curvature of the streamlines. 
Convex curvature decreased both the length and velocity scales of 
turbulent motions, whereas concave curvature caused the opposite. 

The effect of small curvature was found to be much larger than 
one would expect from a linear interpolation between the effects 
of zero and strong curvatures. That is to say, the effect is 
nonlinear, and is strong even at small curvatures. They also ob- 
served that curvature had a relatively larger influence on the 
turbulent shear stress than on the turbulent kinetic energy. An- 
other important conclusion was that convex curvature had a stronger 
effect on the behaviour of the boundary layer than concave curva- 
ture of the same magnitude, (this conclusion is again important 
in the understanding of impeller flows, especially the onset of 
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the wake on the convex suction side of the blades). In addition 
Shivaprasad and Ramaprian provided a way to calculate these ef- 
fects using a correction factor to existing phonoraeno lo gi cal 
me thods . 

A simple description can be given for the mechanism of the 
curvature effect. A flow is considered to be stable* if a fluid 
particle, on being perturbed perpendicularly to the average flow, 
encounters a net restoring force and decelerates. In flows over 
curved surfaces, the centrifugal force is largely balanced by a 
normal pressure gradient. If a fluid particle is displaced away 
from the wall into a region of higher mean velocity, its move- 
ment in this direction is on the average reduced if it flows along 
a convex surface. This is because in its perturbed position the 
centrifugal force acting on the fluid particle will be smaller than 
the mean normal pressure force existing in the flow field. The 
result is a net restoring force. The opposite is true for flow 
along concave walls. Similar considerations also hold for particle 
displacement towards the wall. The effects of rotation on turbu- 
lence are similar, with the difference that the centrifugal ac- 
celeration due to curvature is replaced by the component of the 
Coriolis acceleration normal to the blade surface. 

4.2 Effects of Rotation : As in flow along curved streamlines, 

the influence of rotation is a result of the relation between pres- 
sure forces and inertia forces in the flow. Here the inertia forces 



^Stable here refers to the stability of turbulent fluctuations and 
not the stability of the averaged flow to resist possible separation. 
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are caused by the Coriolis acceleration. In turbomachines the 
Coriolis acceleration has only two components, radial and circum- 
ferential. The circumferential component, which causes the main 
effect of rotation on turbulence since it is almost perpendicular 
to the blade surface, is very small in axial flow machines because 
the radial component of the relative velocity is small. In centri- 
fugal machines, on the other hand, the circumferential Coriolis 
component has a considerable magnitude, and therefore effects of 
rotation on turbulence play a major role. Johnston [84] pointed 
out that there are two basic effects associated with boundary 
layers on rotating surfaces: 1) If a component of the Coriolis ac- 

celeration exists parallel to the solid wall, secondary flows will 
tend to develop in the mean flow field of the boundary layer. This 
effect is present both in axial flow and in radial flow machines 
2) If a component of the Coriolis acceleration exists perpendicular 
to a solid surface, damping or amplifying effects are observed in 
the structure of turbulence. Because of geometrical reasons this 
effect is not very strong in axial machines. In radial machines, 
however, it is of dominant importance. 

In 1972 Johnston, Halleen, and Lezius [85] published experi- 
mental results from a slowly rotating water channel. They con- 
cluded that three different effects can be distinguished in rotat- 
ing boundary layers; 1) The reduction of turbulent wall layer 
streak bursting rate in locally stabilized layers. In a centri- 
fugal impeller these layers are on the suction side of the rotat- 
ing passage. The opposite is true for locally destabilized layers, 
which in a centrifugal impeller are located on the pressure side; 
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2) Production of turbulence is totally suppressed in locally 
stabilized layers; 3) Roll cells are developed on the destabil- 
ized side of the channel. These three mechanisms were correlated 
to the gradient Richardson number and to the rotation number as 
defined in [85]. Later, in 1973, Johnston extended his observa- 
tions to derive a predictive method for this type of boundary 
layer flow [86]. Using the Eddy Reynolds number, Johnston tried 
to predict the rotating boundary layers using a mixing length 
correction that accounts for rotation. Lack of sufficient data 
prevented verification of his results. On the other hand, ap- 
plication of the Eddy Reynolds number criterion to the prediction 
of transition, also observed experimentally, was quite successful. 
Johnston’s results were derived for rotating passages of constant 
cross section. 

To approximate more closely the conditions in centrifugal im- 
pellers, where the passage cross section increases in the flow 
direction, Rothe and Johnston carried out experiments in rotating 
diffusers [87], Their experiments demonstrated that rotation 
strongly enhanced the tendency of a diffuser to stall. This stall, 
however, appeared only on the suction side of the rotating diffuser, 
and was highly steady and two dimensional. The appearance of the 
suction side stall was correlated by Rothe and Johnston to the 
rotation number, and a stall regime map was drawn to enable stall 
prediction in rotating diffusers. It is noted that the appearance 
of stall on the suction side only is fully consistent with previous 
observations of the effect of rotation on the structure of turbulence 
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and is explained by them. In 1976 Rothe and Johnston published 
more data to provide a quantitative account of the onset of flow 
separation in rotating diffusers [88]. The conclusions were simi- 
lar to the conclusions in the previous work [87]. It was, however, 
emphasized in addition that the suppression of turbulent mixing 
and shear stress by the Coriolis acceleration was powerful even 
when fully turbulent upstream wall layers were present, and even 
at small rotation numbers relative to those typical of centrifugal 
impellers. This conclusion is analogous to the observation that 
even very small streamline curvatures have a relatively strong 
effect on turbulent viscosity. 

Very recently Koyama, et al , [89] reported a study confirming 

the observations made thus far. They contributed additional in- 
formation on the quantitative effect of rotation on the skin 
friction coefficient, and on the Monin Oboukhov coefficient, to 
enable more reliable predictions. 

4 . 3 Combined Effects of Surface Curvature and Rotation and 
Their Application : 

So far the effects of streamline curvature without rotation 
and of rotation without streamline curvature were discussed. In 
reality, however, they occur together and must be so treated. The 
effects of the impeller rotation and the curvature of its blades 
on the stability and structure of the internal boundary layers 
must be properly described if the flow pattern in a centrifugal 
impeller is to be predicted correctly. Especially affected by 
curvature and rotation is the onset of boundary layer separation 
and stall in the impeller. The regions of stalled or low velocity 
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flow (wakes) reduce efficiency, reduce operating range, and dis- 
tort the flow delivered at the impeller outlet. 

Johnston and Eide [90] devised a method, using previous ex- 
perimental information, to include the effects of curvature and 
rotation in an existing differential method for turbulent boundary 
layer calculation; namely a slightly modified version [91] of 
STAN 5 [92], STAN 5 was based on Patankar and Spalding’s comput- 

ing scheme [93]. Johnston and Eide modified the turbulent vis- 
cosity model in STAN 5 to include the effects of curvature and 
rotation. They assumed that the mixing length for the boundary 
layer on a flat wall with zero rotation must be multiplied by a 
correction factor F given by 



F ■ 1 - S C Ri c - e a Ri n 
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where (3^ and 3^ are empirical constants for curvature and ro- 
* tational effects y R± are the corresponding Richardson numbers. This 
simple correction factor, which assumes that the effects of curva- 
ture and rotation are linearly combined, was tested on a number 
of boundary layers. Unfor tunately ? all the available cases were 
either curved or rotating and the combination of the two effects 
could not be checked. The agreement between the description 
provided by Eq . (1) and the experiments was good in cases where 

the thin-shear-layer approximation was applicable. 

In 1977 Launder, Priddin and Sharma [94] published a method 
which was applied to spinning cones. It is nevertheless of in- 
terest here. Again Patankar and Spalding’s approach [93] was used 
for the boundary layer calculation in an adapted version of the 
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GENMIX program. They corrected the turbulent viscosity for curva- 
ture and rotation effects using a model of trubulence kinetic 
energy and its local rate of dissipation. The direct effect of 
curvature in the model was limited to a single empirical coeffi- 
cient whose magnitude was directly proportional to a Richardson 
number based on a time scale of the energy-containing eddies. 

They carried out computations for shear flows on curved or spin- 
ning surfaces. They found that the turbulent kinetic energy- 
dissipation rate model did describe the experimental results, pro- 
vided that the transport equation for dissipation was modified to 
include an extra term to account for curvature and swirl effects. 
The success on curved surfaces was better than on spinning sur- 
faces, where the predictions obtained were only marginally super- 
ior to those obtained with a mixing length model [95]. They con- 
cluded that when a mixing length profile can be predicted with 
reasonable certainty, its use would be advantageous for spinning 
flows because of the smaller computing time required to obtain a 
solution. 

The present knowledge on turbulence in rotating passages al- 
lows, at least, qualitative understanding of many of the internal 
passage flow mechanisms. It can explain the occurence of the jet- 
wake flow. It can explain why inviscid pressure side separations 
like the one shown in figure 7 were never detected experimentally. 
It can also provide qualitative justification for observed three 
dimensional flow effects. However, the present knowledge is not 
sufficient to allow accurate quantitative predictions to be made 



46 



of the turbulent viscosity and incorporation of those predictions 
in the 3-D viscous computing techniques now under development. 

Computing techniques, once reliable, would allow for example 
a detailed study to be carried out of the influence of blade sweep- 
back and curvature on stage efficiency and range of operation. 
Presently, the effects of blade sweepback are qualitative infer- 
red from an examination of global overall information such as 
stage performance maps, which correlate blade geometry with ef- 
ficiency and range, but always leave the internal flow details 
in the dark. A detailed analytical study could, for example, ex- 
amine the apparently large and harmful effects of even a small 
amount of convex suction-side blade curvature on the onset of jet- 
wake flow, and might correlate this information with observations 
of the triggering of surge. Also, viscous computations with cor- 
rect viscosity models could be used to create diagrams like the one 
in Fig. 1, and give a better understanding of the processes caus- 
ing jet-wake flow. 

5. Inner-Outer Solution Combinations (Patching Techniques) 

A standard approach to solving for a complete flow field, 
is fir"st to divide the flow field into regions. In each region, 
a different simplifying assumption is made in order to allow an 
easy solution to be obtained. The regions are connected by com- 
mon boundary conditions on common boundaries. The classical ex- 
ample of this approach is boundary layer theory. In centrifugal 
impellers a similar approach was used by a number of investigators 
in a very simplified way to attack the problem of the jet-wake 



flow* ** The problem is of course complicated by the fact that the 
region boundaries are not known a-priori. Consequently, the tech- 
niques are either iterative or depend strongly on empirical re- 
sults. At present, the information quoted above on turbulent vis- 
cosity variation in curved rotating passages is not directly in- 
corporated in viscous flow field calculations, and some methods 
avoid viscosity altogether. 

In 1973 John Moore [96] published a method of calculating the 
jet-wake flow in a rotating radial flow passage. Moore divided 
the flow field into four parts; 1) Potential flow; 2) Top and 
bottom wall boundary layers; 3) Corner flows, and 4) Side wall 
boundary layers. Moore calculated the potential core flow ac- 
counting for the displacement and acceleration of the potential 
flow by the wake on the suction side wall. The development of the 
cross flows on the top and bottom walls was computed on the basis 
of Moses 1 strip integral method [97] modified to include the three 
dimensional and rotational terms for steady incompressible flow 
in rotating passages. Thus on top and bottom the wall boundary 
layers as well as radial and tangential flows were predicted. 

The corner flows were modeled by assuming that continuity holds 
over the corner cross section, and that the fluid turning through 
the corner maintains its momentum in the radial direction. The 
side wall boundary layers were computed with an analysis which 
included a Coriolis 1 acceleration term associated with velocities 



*There are also a number of works published on boundary layers and 
secondary flows in impellers. Though relevant, these studies are 

not included in order to keep this review within reasonable limits. 
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normal to the wall and a Coriolis term affecting the pressure 
gradient normal to the wall. The model gave a momentum integral 
equation similar to the one derived by Moon [98]. However, it 
was applied with irrotational potential flow at the edge of the 
boundary layer, and retained the terms which arose due to the 
two Coriolis effects. Moore calculated the wake assuming that it 
started at a point where the shape factor was 3.0. Subsequent 
flow was considered to consist of similarity profiles, including 
a cross flow entrainment into the wake from the corners. These 
four partial flows were combined to yield a good agreement between 
calculations and experiment (Fig. 17). 

In 1975 Sturge and Cumpsty developed a two dimensional method 
to predict the incompressible jet-wake flow [99]. For simplicity, 
they avoided the direct modelling of viscous effects, which led to 
major drawbacks in their method. The method was later extended 
to compressible flows [100]. The two dimensional blade to blade 
passage was divided into two regions; 1) a jet region where the 
flow was assumed to be inviscid and irrotational; 2) a wake region 
where there was no flow. The boundary between the two regions 
was treated as a free shear layer without mixing between the two 
regions, similar to an air-water interface. The position of the 
shear layer was characterised by the condition that there was no 
change in static pressure across it, and that in the wake the pres- 
sure was governed only by hydrostatic effects due to the centri- 
fugal acceleration field. The flow in the jet was analysed using 
irrotat ionali ty and continuity and defining a stream function, 
and was solved numerically with the shear layer between the jet 
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and the wake as one of the boundaries. This required an iterative 
process which was shortened by a variable transformation which left 
the stream function as one of the independent variables. The 
resulting second order partial differential equation was solved 
using a finite difference technique. Besides the fact that the 
method was restricted to irrotational flow, a much more disturbing 
feature was that the separation point on the suction side could 
not be predicted, but had to be known a-priori or assumed. (This 
was the penalty for not modelling the viscous effects). Though 
a separation criterion based on a velocity ratio between separa- 
tion and the highest velocity attained on the suction side was 
suggested by Sturge and Cumpsty the method was still not predict- 
ive unless the calculation of a separation point on the suction 
side was incorporated. Another drawback of the technique was that 
no flow was assumed in the wake. This is not in accordance with 
Eckardt ’ s observations [12]. Eckardt measured significant veloc- 
ities inside the wake. The magnitude of the velocities was about 
30% of the jet velocities. No comparison to experiments were given 
in Sturge and Cumpsty ’s work. 

Howard and Osborne [101] developed a further simplified method 
to predict the jet-wake flow. The method was based on Katsanis ' 
streamline curvature inviscid method on the H-S stream surface [21] 
combined with an extremely simplified calculation on the B-B sur- 
face [102]. In contrast to Sturge and Cumpsty f s method the analy- 
sis did not determine an exit flow angle and required a separate 
correlation for this purpose. The wake was treated, as in the 
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previous method, as an additional blockage in the passage. The 
inviscid calculation was carried out in the jet only with the mid 
passage H-S surface, moved towards the pressure side by one-half 
the wake width. Subsequently, the B-B analysis was carried out in 
the region of the jet only. The continuity equation at each 
quas i-or thogonal was adjusted to account for growth of the wake 
width. The wake width and the jet flow angle , required in the cal- 
culation, were not predicted by the technique but had to be estab- 
lished either by empirical correlations or from measured data. 

The ability of the method to be really predictive depended on the 
accuracy and the generality of this empirical information. Al- 
though some three dimensional aspects of the jet-wake flow were 
included, they were possibly oversimplified by the assumption 
that the wake was evenly distributed along the suction side of the 
passage from shroud to hub. This is not consistent with Eckardt 1 s 
observations [12] . 

The problem of calculating the jet-wake flow properly is form- 
idable. Only three attempts have been reported so far, and none 
was based on a partial differential approach. Two (Sturge and 
Cumpsty, and Howard and Osborne) avoided viscosity altogether by 
replacing the wake with a region without flow. Moore, in a more 
detailed model, applied an integral boundary layer technique. 

The full task is still not accomplished and further attempts are 
anticipated . 

6. Fully Viscous Solutions 

A large volume of work on the solution of viscous flows was 
published in the past few years and a complete review will not be 



51 



attempted. Of direct interest here are solutions for flows which 
are both rotating and turbulent. A number of relevant works will 
be chronologically reviewed. 

In 1975 Walitt, Harp and Liu published a technique [103], to 
predict the viscous flow field on the B-B stream surface. It was 
the first phase of a technique developed to solve the steady three 
dimensional compressible Na v i er - S to ke s equations in rotating curvi- 
linear coordinates. Any turbulence model could be incorporated in 
the method which was applied (in. Phase I) to laminar flow and to 
turbulent flow described by a simple mixing length correlation not 
incorporating rotation and curvature effects. In 1978 the second 
phase of the work, solving the flow on cross sectional surfaces 
perpendicular to the direction of the primary flow, was published 
[104]. The full technique solved the equation of motion by itera- 
tions. The solution started from a zeroth iteration which was a 
3-D inviscid solution; Katsanis ’ method [21] was used. The solution 
marched alternately from hub to shroud (solving on B-B surfaces) 
and from inlet to outlet (solving on cross-sectional surfaces). 

The cross sectional iteration employed the B-B flow field as the 
previous iterate. The method was an explicit two dimensional time 
marching scheme extended to three dimensional flows using the 
equivalence principle of Hayes. In its original derivation, the 
equivalence principle states that for slender bodies at hypersonic 
speeds the three dimensional steady equations of motion for in- 
viscid flow reduce identically to unsteady equations in two space 
dimensions. An extension of the principle was used in forming an 
analytical model which calculated the viscous cross flow with a 
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known flow from a previous iterate in the marching direction. 

The coordinate in the marching direction was made proportional 
to a time like variable by forming the product with a character- 
is tic velocity . 

The method was used to predict the flows in a radial exit, and 
in a backswept compressor impeller. In the radial exit impeller a 
separation zone was predicted on the suction side of the passage 
close to the outlet. The flow in the backswept impeller was at- 
tached. No conclusions can be drawn with regard to the effective- 
ness of sweepback in reducing or eliminating the jet-wake flow 
because the radial exit impeller was calculated with a laminar 
viscosity model, and the backswept impeller was calculated with a 
turbulent mixing length model for steady flows along straight 
surfaces. No comparison to experiments was given, but in [104] a 
qualitatively correct jet-wake pattern was predicted. 

In 1976 Briley, Kreskovsky and McDonald published a method 
developed for stationary axial trubomachinery passages [105]. It 
is included here because it can in principle be extended to ro- 
tating centrifugal passages. The method was based on an approxi- 
mate formulation which was solved in a forward marching technique 
in a primary spacial direction. The complete flow was obtained 
in a sequence of two dimensional calculations, resulting in a 
substantial saving of computer time and storage over that which 
would be required for the solution of the full, elliptic Navier- 
Stokes equations. Essential to the method was the derivation of a 
particular coordinate system. One coordinate direction was 



identified with the primary flow direction along which the solution 
was to march, while the two remaining coordinate directions could 
be associated with secondary flows. The full elliptic flow model 
was simplified by a number of assumptions. For high Reynolds 
numbers, viscous effects were negligible except in the boundary 
layers on all four walls. The method was an extension of 3-D 
boundary layer methods, but the approximate equations were also 
used in the inviscid core flow. For entirely supersonic flows 
the equations, together with initial and boundary conditions, could 
be solved by forward marching in the primary direction without any 
assumptions about the pressure field. For subsonic, however, where 
the inviscid core was elliptic, downstream boundary conditions 
were required to enable a solution. This problem was overcome by 
imposing on the solution, an entirely inviscid pressure field which 
was calculated a priori, with proper elliptic boundary conditions. 

The pressure field then had to be corrected as part of the solution 
procedure. The viscous pressure drop in the primary direction was 
treated separately as a function of the distance in this direction 
only, and was computed with the forward marching integration using 
integral mass conservations in that direction. No assumptions were 
made concerning the pressure gradients in the secondary flow direc- 
tions. The turbulence model incorporated in the method was based 
on the solution of conservation equations for turbulent kinetic 
energy and dissipation, and as the technique was developed for 
stationary passages, rotation effects were not included. Curvature 
effects were also neglected. 



Another viscous computation technique for stationary turbo- 
machinery ducts was reported by Dodge [106]. Again the model was 
not based on the full Navier- S to kes equations in order to save 
computer time and storage space. The method was based on separat- 
ing the momentum and continuity equations into coupled elliptic and 
parabolic parts which were then solved iteratively. Dodge sepa- 
rated the velocity into a viscous and a potential component. Sub- 
stitution of the decomposed velocity into the equation of motion 
yielded a simpler equation of motion from which the pressure gradi- 
ent was eliminated using the velocity potential. The velocity 
potential in turn was expressed using the continuity equation. 

The simplified equation of motion was parabolic with a primary 
direction in the flow direction, and was solved using a marching 
technique similar to those used in boundary layer theories. The 
equation for the potential was elliptic and was solved with a 
relaxation technique. Since both equations contained coupling 
terms, iteration was required between the two. Dodge developed and 
used his technique for stationary passages; consequently his turbu- 
lence model did not include effects of rotation. He also did not 
correct the turbulence model for curvature effects despite the 
presence of curvature in the passages which were analysed. The 
turbulent viscosity was calculated using a mixing length concept 
near the walls (including Van Driest ' s wall shear correction and 
Cebeci-Smith pressure correction) and a constant free stream vis- 
cosity away from the walls. A nearly constant viscosity was used 
in wakes. Dodge ? s technique failed in the presence of reversed 
flow, the forward marching differences then becoming unstable. 
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By replacing the variation of the dependent variable in the normal 
direction by a spectral series [107] this problem was overcome. A 
number of comparisons with the experiments of Schubauer and Klebanoff 
revealed good agreement. 

In 19 7 7 a paper [108] was published by Majudmar , Pratap and 
Spalding which described a method based on Patankar and Spalding f s 
boundary layer calculation [109]. The method of Patankar and Spalding 
for parabolic flows is based on the uncoupling of the pressure; 
using in the momentum equation for the primary direction (the radial 
direction) a pressure value which is averaged over the cross section. 
In the method of Majudmar et al the model reported in [109] was 
modified through addition of a Coriolis term in the circumferential 
momentum equation and both a centrifugal and a Coriolis term in 
the radial momentum equation. The turbulence model used was not 
corrected for rotational and curvature effects; it was based on 
Launder ! s and Spalding's kinetic ener gy- d i s s ip a t i on concept [110]. 

The model was solved in a single forward marching sweep from inlet 
to outlet in the following steps: 

1) The pressure distribution in a cross section perpendicular 
to the primary direction, and the average value of the 
pressure in that cross section were assigned preliminary 
gues sed values . 

2) The three momentum equations at a downstream position 
were solved, using flow properties at an upstream posi- 
tion, to yield a first approximation to the velocity 
distribution . 
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The mean cross sectional pressure was corrected to give 
new radial velocity components which satisfied continu- 
ity. This was done using a linearized radial momentum 
equation. 

4) A corrected pressure field on the cross section surface 
was computed using an elliptic equation derived from 
continuity, and two linearized momentum equations in 
that surface. The cross sectional velocities were cor- 
rected according to this pressure field. 

5) A new turbulent viscosity distribution was evaluated from 

a solution of the kinetic energy and dissipation equations. 
Figure 18 illustrates the good agreement obtained between 
this technique and Moore’s experiments [10]. Majumdar, Pratap and 
Spalding concluded, however, that their parabolic technique was 
only accurate for low values of QB/w (where £2 is the angular 
velocity, B is the circumferential width of passage, and w is the 
bulk mean velocity). For high values of ftB/w, partially para- 
bolic procedures [111, 112] in which the pressure uncoupling was 

abandoned yielded more accurate results. Both the partially para- 
bolic technique [111, 112] and the parabolic technique [108] were 

based on identical mathematical and turbulence models. The dif- 
ference was in the method of solution. While in the parabolic 
technique a single marching sweep was carried out from inlet to 
outlet, the partially parabolic technique was necessarily itera- 
tive. A number of marching sweeps, each similar to the single 
sweep of the parabolic method, were required. In each sweep, or 
iteration, the pressure field was corrected until convergence was 
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achieved. Thus through the pressure field, downstream effects 
were transmitted into the flow field and the solution was given a 
radially, or streamwise, elliptic character. Results for secondary 
flows were given in [112] and compared to measurements of Wagner 
and Velkoff [113]. Agreement was good in some cases, but signifi- 
cant disagreements were found in others. 

Very recently Bosman, Chan and Hatton published a finite dif- 
ference technique which was simple and probably not expensive to 
run [114]. The model was, however, subjected to a number of restrict- 
ions which limited its application to two dimensional flows in 
radial surfaces which were absolutely plane and perpendicular to 
the axis of rotation. The flow had to be incompressible and the 
viscosity constant in the entire flow field. Further, the tech- 
nique utilized a nond imens iona 1 slip parameter to represent the 
blade surface shear stresses. The slip parameter was assumed 
to be constant along the blades from inlet to outlet, and reliable 
information was not given on its evaluation. Using the incompres- 
sible equations of continuity and momentum and a stream function, 
the flow was formulated with two coupled partial differential el- 
liptic equations. This was the distinctive feature of this tech- 
nique in contrast to the parabolic, or partially parabolic nature 
of the forward marching methods. The unknown functions to be 
evaluated were the stream function and the relative vorticity. 

The equations were solved with a relaxation technique with finite 
difference approximations. Central differences were used for dif- 
fusion terms and up-wind differences for convection terms. The 
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two equations were solved iteratively until convergence was achieved* 
The results were compared to measurements carried out in a water 
[115, 116] and to results of an inviscid calculation based on 

Wu. At the design point of the impeller, agreement between the 
experiments, the inviscid calculation, and the viscous calculation 
was not bad. This is not surprising in view of the fact that the 
impeller was two dimensional and strongly backswept. In such a 
configuration, viscous effects are not very strong (see Fig. 1) and 
inviscid calculations should give fairly accurate results. The 
situation at off design points, with positive or negative inlet 
incidence angles, was not much different downstream and close to 
the impeller exit. At upstream positions, closer to the inlet, 
the agreement between experiments and viscous results was better 
than agreement between experiments and inviscid results. This 
again is not surprising in view of the fact that inviscid theories 
should not be used at all at off design conditions since they can 
not represent properly the viscous effects associated with the 
inlet shock. The technique of Bosman et al [114], when further 
developed to remove some of its restrictions could become a useful 
design aid . 

Also in 1978 Goulas and Baker [117] published a method for 
calculating the viscous flow on either the H-S or on the B-B 
stream surfaces. The method was not three dimensional thus far, 
but the author f s s t a t e d in t en t ion was to couple the two solutions 
in an iterative sequence to obtain the three dimensional flow field. 
The published method was based on the usual conservation equations 
written for a stream surface. Introducing a stream function, the 



two equations of motion for the stream surface were combined into 
a Poisson like equation, V 2 \p = q(r.cj).Z), in which the function 
q was not known. The function q contained all the differences 
between incompressible irrotational flow and real viscous compres- 
sible flow under consideration. The solution by finite differences 
on a stream surface was carried out in the following steps: 

1) A stream function distribution was assumed from which the 
initial two components of the velocity on the surface 
were calculated. The third velocity component was cal- 
culated using the geometry of the stream surface. 

2) With this information the turbulent viscosity was computed 
using Launder ? s and Spalding’s kinetic energy-dissipation 
model [110]. A set of linear equations was formed and 
solved simultaneously for each quas i-or thogonal line. 
Starting from the inlet where the kinetic energy and 
dissipation were known as initial values, the solution 
proceeded downstream to the exit. 

3) The streamwise momentum equation was integrated along 
streamlines, using the calculated turbulent viscosity, 
to yield an entropy field. 

4) With velocities, entropy and density known, the function 
q was calculated, and Poisson equation for the stream 
function was then solved. 

The process was iteratively repeated until convergence of the 
stream function was achieved. 
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Results were compared with the inviscid solution of the flow 
field in a turbocharger compressor and viscous effects were qual- 
itatively clearly visible [117]. No comparison with experiments 
was given. 

7. Comparison of Computation Times 

The solution techniques mentioned in this review require 
numerical methods. A comparison of their drawbacks and merits is 
not complete without consideration of their demands for computer 
time and computer storage. For most numerical techniques the 
storage requirements are not published, though the particular com- 
puter model used can sometimes indicate what might be required. 

The following table gives an approximate indication of computation 
times for some of the techniques mentioned in this review. An 
exact accounting, however, is not possible as some authors mention 
virtual execution time, and others indicate a total time which in- 
cludes compilation of the program. Some authors even neglect to 
define the time that they quote. Furthermore, the time stated 
depends strongly on the problem solved and the accuracy obtained 
in the results. The accuracy, in turn, is affected by mesh size 
and convergence criterions adopted. The values given in the 
following table should therefore, be treated as an order of mag- 
nitude indication only. Table I classifies the techniques accord- 
ing to the following abbreviations: 

SC - Streamline Curvature 
FD - Finite Differences 
FE - Finite Elements 
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Table 1. Comparison of Approximate Computation Times 



Author 


Ref. No. 


Computer Used 


Katsanis 


21 


IBM 360 


Marsh 


27 


KDF 9 


Katsanis /McNally 


29 


IBM 360 


Katsanis/McNally 


30 


IBM 360 


Ad 1 e r / Kr im er ma n 


31 


IBM 370 


Frost 


34 


SDS 90/300 


Davis/Millar 


35 


? 


Davis /Millar 


35 


? 


Wilkinson 


43 


IBM 360 


Katsanis/McNally 


49 


IBM 360 


Adler/ Kr im erman 


54 


IBM 370 


Katsanis 


56 


IBM 360 


Senoo/Nakase 


68 


Facom 230-60 


Katsanis 


66 


IBM 360 


Senoo/Nakase 


67 


Facom 230-60 


Bosman/Shaarawi 


69 


CDC 7600 


Novak/Hearsey 


70 


CDC 6600 


Hirsch/Warzee 


73 


CDC 6500 


Kr imerman/ Adler 


74 


IBM 370 


Wallit/Har p/Liu 


103 


CDC 6400 


Dodge 


106 


CDC 6400 


Dodge/Lieber 


109 


CDC 6400 


Dodge/Lieber 


109 


CDC 7600 


Madjumar/Pratap 


108 


CDC 6600 


Mad j umar/ Spalding 


111 


CDC 6600 



Time Flow Techniqu 

(Mins) 



10 


H-S 


inviscid 




SC 


10-15 


H-S 


inviscid 




FD 


30 


H-S 


inviscid 




FD 


3-15 


H-S 


inviscid 




FD 


1.5 


H-S 


inviscid 




FE 


10-15 


H-S 


inviscid 




SC 


1 


H-S 


inviscid 




SC 


.5 


H-S 


inviscid 




FD 


.5 


B-B 


inviscid 




SC 


30 


B-B 


inviscid 




FD 


2 


B-B 


inviscid 




FE 


10 


B-B 


transonic/inviscid SC/F 


10(sec) 


B-B 


inviscid 


(cpu) 


FD 


2 


3-D 


inviscid 




SC/F 


16 


3-D 


inviscid 




SC/F 


3.5 


3-D 


inviscid 




FD 


1.5 


3-D 


inviscid 




SC 


1 


3-D 


inviscid 


(cpu) 


FE 


5 


3-D 


inviscid 




FE 


3000 


3-D 


viscous 






60 


3-D 


viscous 






7 


3-D 


viscous 






.5 


3-D 


viscous 






2.5 


3-D 


viscous 






2.5 


3-D 


viscous 
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Figure 1: Possible flow regime classification in centrifugal impellers 

(the location and width of the border zone are assumed and 
have yet to be verified) . 
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axial turbine taken from Ref. 34. 
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axial total total axial total total 

VELOCITY PRESSURE TEMPERATURE VELOCITY PRESSURE TEMPERATURE 

m/s kPo °K m/s kPa °K 
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Figure 4: Comparison between experiments, a H-S matrix method and a H-S 

streamline curvature technique according to Davis and Millar 
[35]. The example is an axial compressor taken from Ref. 35. 
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Figure 5: Comparison between experiments and Hirsch andWarzee's H-S 

finite elements calculations [33]. The example is an axial 
compressor taken from Ref. 33. 
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STATIC PRESSURE COEFFICIENT MACH NUMBER 




o) STATOR 




NORMALIZED DISTANCE ALONG PASSAGE 
b) TURBINE BLADE 

Figure 6: Comparison between experiments and Wilkinson's B-B 

streamline curvature calculations [43]. The examples 
are an axial compressor and an axial turbine taken 
from Ref. 43. 
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Figure 7: Streamlines inside the rotating impeller as 

calculated by Stanitz [45]. The predicted 
inviscid pressure side separation is well 
visible on the right side. Taken from Ref. 45. 
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Figure 8: Comparison between experiments, Stanitz f s B-B calculations and 

Adler and Krimerman’s B-B calculations [54]. The example is 
a centrifugal impeller taken from Ref. 54. 
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[55]. The example is an axial turbine stator taken from Ref. 55. 
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Figure 10: Comparison between transonic experiments 

and Kurzrock and Novick's B-B calculations 
[62], The example is an axial compressor 
cascade taken from Ref. 62. 
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Figure 11: Comparison between experiments and Senoo and Nakase's 

quasi 3-D calculations 167]. The example is a 
centrifugal impeller taken from Ref. 18. 
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Figure 12: The difference between the traces of 2-D surfaces (dashed 

lines) and 3-D stream surface traces, calculated by 
Krimerman and Adler [74]. The example is a centrifugal 
impeller taken from Ref. 74. (the section 6 to 9 shown 
are located in the inducer and beyond, in the region where 
the flew undergoes a double-change of direction both in 
the tangential direction and from axial to radial) . 
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PRESSURE COEFFICIENT 




Figure 14: Comparison between experiments and 3-D inviscid 

calculations by Krimerman and Adler [74]. The 
example is a centrifugal radial exit impeller 
taken from Ref. 74. 
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Figure 15: Comparison between Eckardt's measurements at 

the passage center [12] and 3-D calculations 
by Krimerman and Adler [74]. ( a-immediately 

after inlet; b,c-about mid passage location; 
d-close to exit) 
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Figure 16: Comparison between Eckardt's measurements near 

the Shroud and 3-D calculations by Krimerman 
and Adler [74]. (b,c-about mid passage 
location) 
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Figure 18: Comparison between the experiments of Moore [10] and Majumdar f s et. al., calculations [108]. 
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